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I.  INTRODUCTION 

A.  MOTIVATION 

Accurate  modeling  of  progressive  flooding  is  essential  in  assessing  the  survivability 
features  of  a  given  ship.  A  realistic  flooding  and  counter  flooding  scenario  is  a  very 
complex  dynamic  system  including  complicated  vessel  motions,  time-varying  external 
excitation,  changing  hull  mass  and  hydrodynamic  properties,  progressive  flooding 
through  damage  openings,  internal  sloshing  of  flooded  water,  and  de-watering  control 
actions.  Although,  due  to  the  complexity  of  the  problem,  the  traditional  method  of 
approach  is  through  large  amplitude  static  analyses,  the  need  for  a  more  accurate  analysis 
has  been  recognized.  In  order  to  bridge  this  gap,  recently  there  has  been  a  large  amount  of 
research  focusing  on  dynamic  stability  prediction  taking  into  consideration  water-on-deck 
and  large  amplitude  time-dependent  motions  and  excitation  [Ref.  1].  Equally  important, 
however,  is  the  development  of  a  design  methodology  for  incorporating  elements  of 
man/machine  interface  into  the  problem;  in  other  words  effects  of  decision-making  and 
control  strategies.  This  presents  an  additional  level  of  difficulty  since  combinations  of 
discrete  events  and  continuous  differential  equations  describe  such  systems.  The  long 
term  goal  that  is  initiated  with  this  thesis  is  to  study  the  development  of  an  improved 
damage  stability  analysis  and  design  tool,  which  will  address  the  effects  of  control 
decisions  and  vessel  dynamics  in  a  progressive  flooding  scenario. 


B.  BACKGROUND 

When  a  ship  suffers  underwater  damage,  there  is  a  rapid  influx  of  flooding  water, 
usually  followed  by  a  period  of  slower  progressive  flooding.  Progressive  flooding  is 
defined  as  flooding  which  extends  beyond  the  compartments  that  are  open  directly  to  the 
sea.  This  results  in  flooding  of  compartments  whose  hull  boundaries,  but  not  interior 
bulkheads,  are  still  intact.  Existing  computer  models  at  the  Naval  Postgraduate  School 
utilize  the  FORTRAN  computer  language  and  formatted  input  files  to  model  progressive 
flooding  [Ref.  2].  Any  change  in  model  design  conditions,  flooding  damage,  or  damage 
control  abilities  requires  a  completely  new  FORTRAN  input  file.  The  purpose  of  this 
thesis  is  to  include  the  capabilities  of  the  original  computer  progressive  flooding  model  in 
a  more  modular  and  expandable  format,  using  the  user-friendly  MATLAB  and 
SIMULINK  programs. 

C.  DEFINING  THE  MODEL 

1.  Hole  Model 

Holes  can  be  modeled  as  a  short-tube  orifice  as  shown  in  figure  1  with  a  diameter 
(hole  size)  much  larger  than  its  tube  length  (ship's  skin  thickness).  The  flow  discharge 
coefficient  for  this  case  is  0.816  [Ref.3:  pp.42-43]. 

2.  Hull  Model 

The  required  hull  characteristics  are  determined  using  the  same  table  of  offsets 
used  by  Dawson  (Hull.dat)  [Ref.  2:  pp.  127-143]  and  MATLAB  programs  arvol .  m  and 
drafts  .  m  to  calculate  the  locations  of  center  of  buoyancy  and  center  of  floatation, 
waterplane  area,  and  waterplane  moments  used  in  the  program  equations.  Both  Matlab 


programs  are  included  in  the  appendix.  For  the  FFG-7  hull  model  used  here  the 
following  values  were  calculated: 

Table  1 :  Hull  Model  Characteristics 
Center  of  Buoyancy 

Xbar 209.9  ft  aft  of  FP 

Ybar 0.0  ft  from  CL 

Zbar 10.52  ft  from  baseline 

Center  of  Floatation 

Xbar 227.4  ft  aft  of  FP 

Ybar 0.0  ft  from  CL 

Waterplane  area 14030  ft 
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Longitudinal  Moment  of  Inertia 1.5883  x  10  9  ft  4 

Transverse  Moment  of  Inertia 1.8148   x  10  6  ft  4 

3.  Equilibrium  Model 

The  added  weight  method  [Ref.4:  pp. 76-79]  is  used  within  the  program  to 
determine  the  equilibrium  condition  (displacement,  heel,  and  trim)  of  the  hull  model  for 
any  state  of  flooding.  This  method  assumes  hull  characteristics  are  constant  and  is 
accurate  for  "small  angles".  The  added  weight  method  treats  flooding  water  as  added 
weight  and  treats  the  hull  as  if  it  were  intact.  The  effect  of  the  added  weight  is  divided 
into  a  weight,  which  causes  the  ship  to  immerse  evenly  (parallel  sinkage)  and  moments 
which  cause  the  ship  to  list  and  trim  about  its  center  of  floatation.  The  small  angle 
assumption  assumes  that  the  ship  lists  and  trims  about  a  fixed  center  of  floatation.  In 
reality,  since  the  hull  is  not  cylindrical,  the  location  of  the  center  of  floatation  will  change 
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with  changing  list  and  trim.  This  method  allows  much  simplified  calculation  of  hull 
characteristics  compared  with  the  lost  buoyancy  method.  Since  the  weight  added  is  liquid 
water  which  is  free  to  move  as  the  ship  moves  (and  thus  does  not  contribute  to  pitch  and 
roll  resistance),  a  free  surface  correction  is  added  to  the  longitudinal  and  transverse 
moments  of  inertia  to  account  for  the  reduced  resistance  of  the  hull  to  pitch  and  roll. 

4.  Compartment  Model 

The  ship  compartments  are  modeled  as  rectangular  boxes.  Extension  to  more 
realistic  shapes  is  one  of  the  recommendations  for  further  expansions  to  the  capabilities 
of  the  program. 

D.  SIMULINK 

The  SIMULINK  program  provides  a  graphical  representation  of  mathematical 
relationships  and  equations  to  enable  the  user  to  determine  the  response  of  system 
odels.  Graphical  symbols  and  connections  are  used  instead  of  numbers  and 
mathematical  symbols,  allowing  more  intuitive  use  and  understanding  of  complex 
mathematical  relationships  between  system  components.  The  SIMULINK  program  also 
allows  easy  modification  of  input  parameters  and  import  and  export  of  data  from 
MATLAB  programs.  Lastly,  SIMULINK  allows  the  user  to  observe  system  or  individual 
component  behavior  during  the  simulation  through  the  use  of  output  graphs. 

E.  OBJECTIVES 

The  objective  of  this  thesis  is  to  develop  a  progressive  flooding  model  to  utilize  the 
superior  user  interface  and  data  exchange  capabilities  of  the  SIMULINK  program.  The 


program  must  be  able  to  accurately  calculate  flooding  rates  and  the  subsequent  response 
of  the  ship.  The  program  must  also  be  user-friendly,  easily  modified  and  have  the 
potential  for  further  expansion  by  subsequent  users.  To  achieve  these  goals  the  following 
programs  were  written  by  the  author: 
1.  Programs 

a.  Pflood.m 

Pflood.m  is  the  main  program  created  by  the  author.  This  program  calculates  the 
flooding  and  ship  response  caused  by  a  hole  of  user-chosen  size  and  location  in  the 
ship's  hull.  The  flooding  rate,  flooded  volume,  heel  and  trim  are  calculated  and 
updated  each  second.  In  addition,  the  user  can  specify  the  number,  time  of 
activation  and  size  of  pumps  used  for  dewatering  and  the  duration  of  their  use. 

b.  Arvol.m 

Arvol.m  is  a  MATLAB  program  written  by  the  author  to  determine  the  center  of 
buoyancy  of  the  hull  model  defined  by  the  Hull.dat  data  file  using  Simpson's  rule 
integration  techniques. 

c.  Drafts.m 

Drafts. m  is  a  MATLAB  program  written  by  the  author  to  determine  the  waterplane 
area,  center  of  floatation,  and  second  moment  of  inertia  of  the  hull  waterplane  for 
a  given  draft. 
2.  Validation 

Program  performance  was  validated  by  comparing  the  program  output  to  expected 
ship  behavior  for  a  given  flooding  condition. 


II.  SIMULINK  PROGRAM 

The  schematic  representation  of  Pflood.m  is  shown  in  Figure  1 .  For  illustrative  purposes, 
Pflood.m  is  shown  in  its  "expanded"  form.  By  leaving  Pflood.m  in  expanded  form,  any 
function  or  variable  can  be  observed  or  changed  by  the  user.  For  ordinary  use,  the  various 
subroutines  would  be  grouped  into  subsystem  blocks  with  only  the  user-defined  variables 
and  the  outputs  represented  by  discrete  function  blocks.  Placing  the  program  in 
subsystem  block  form  has  the  added  advantage  of  greatly  enhancing  program  readability 
and  user-friendliness.  Specific  details  of  the  individual  block  diagram  elements  are 
included  in  Chapter  HI. 
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Figure  1 


Pflood.m  Simulink  Program 


III.  PROGRAM  AND  SUBROUTINES 

A.  FLOW  RATE  THROUGH  HOLE(S) 

The  flow  through  the  damage,  either  primary  or  secondary,  can  be  modeled  as 
turbulent  flow  through  an  orifice  of  negligible  length,  or  a  hole.  Flow  rate  through  a  hole 
is  a  function  of  the  hole  size  (cross  sectional  area),  hole  shape,  and  pressure  difference 
across  the  hole.  The  hole  size  and  shape  are  related  to  the  discharge  coefficient,  Cd-  A  Cd 
of  0.816  is  used  in  this  application  in  order  to  remain  consistent  with  Dawson's  previous 
computer  flooding  model.  This  number  is  also  consistent  with  theoretical  predictions  of 
flow  through  orifices  as  explained  in  the  previous  chapter  [Ref.  3:  pp.  42-43].  The 
pressure  difference  is  determined  by  the  difference  in  height  of  the  liquid  on  either  side  of 
the  hole.  The  equation  is 


Q  =  Cd  x  Ax^2g  x(hi-  hi) 

where  Q  is  the  volumetric  flow  rate,  A  is  the  cross  sectional  area,  Cd  is  the  discharge 
coefficient,  g  is  the  acceleration  due  to  gravity,  and  h2-hi  is  the  liquid  height  difference 
across  the  hole.  This  height  difference  is  a  function  of  the  level  of  the  flooding  water 
inside  the  flooded  tanks,  as  well  as  the  draft  and  trim  of  the  ship.  The  latter  quantities 
will  have  to  be  calculated  using  the  added  weight  method  as  explained  in  the  previous 
chapter.  The  Simulink  representation  of  the  above  equation  is  shown  in  Figure  2. 
The  program  allows  continuous  adjustment  of  the  values  of  hi  and  h2  to  account  for 
increased  immersion  of  the  exterior  hole  and  increased  compartment  water  height. 
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Figure  2:  SLMULINK  Hole  Flow  Symbology 


B.  NET  FLOW  INTO  A  COMPARTMENT 

Flooding  water  removed  from  the  compartment  by  pump  activation  or  by  flow 
into  adjacent  compartments  is  subtracted  from  the  water  flow  rate  into  the  compartment. 
Two  way  flow  due  to  water  removal  from  one  of  two  connected  flooded  compartments  is 
also  possible,  and  is  accounted  for  by  taking  the  square  root  of  absolute  value  of  the 
liquid  height  difference  and  then  multiplying  the  result  by  the  sign  of  the  height 
difference.  This  method  avoids  program  failure  due  to  trying  to  take  the  square  root  of  a 
negative  number.  This  program  also  continuously  updates  the  values  of  hi  and  h2  to 
account  for  increased  immersion  of  the  exterior  hole  and  increased  compartment  water 
height. 

C.  HEIGHT  OF  WATER  (HOW) 

Once  the  net  volumetric  flow  rate  of  water  is  known,  integrating  the  flow  rate 
results  in  the  net  volume  of  water  added  to  the  compartment.  This  volume,  divided  by  the 
horizontal  area  of  the  compartment,  and  divided  again  by  the  permeability  of  the 
compartment,  gives  the  height  of  water  in  the  compartment.  A  saturation  function  ensures 
that  the  height  of  water  in  the  compartment  is  always  between  zero  and  the  maximum 
height  of  the  compartment.  Figure  3  shows  the  Simulink  height  of  water  function. 

D.  MOMENTS  DUE  TO  ADDED  WATER 

The  added  weight  method  treats  the  added  weight  as  a  combination  of  a  weight 
which  causes  the  ship  to  sink  evenly  (parallel  sinkage)  and  transverse  and  longitudinal 
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Figure  3:  SIMULINK  Height  of  Water  Terminology 
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moments  which  cause  the  ship  to  heel  and  trim  respectively.  Multiplying  the  height  of 

water  (HOW)  by  the  product  of  the  weight  density  of  the  liquid  ( pg  ),  the  compartment 

permeability  (per),  and  the  compartment  horizontal  area  (A)  determines  the  weight  added 

to  the  compartment.  Multiplying  the  added  weight  by  the  moment  arm  (Ztrans/iong)  from 

the  ship's  center  of  floatation  to  the  center  of  gravity  of  the  added  weight  gives  the 

moment  about  the  center  of  floatation  in  the  longitudinal  and  transverse  directions.  In 

equation  form: 

Moment  =  {pgA  x  HOW  x  per)  x  Ztram  /  long 

The  Simulink  representation  of  the  above  equation  is  shown  in  Figure  4. 

E.  KG  (VERTICAL  POSITION  OF  CENTER  OF  GRAVITY) 

KG  (Vertical  Position  of  Center  of  Gravity)  is  determined  by  multiplying  the 
original  weight  of  the  ship  by  the  original  KG,  adding  the  product  of  the  weight  of  the 
added  water  and  the  height  of  the  added  water  center  of  gravity  above  the  keel,  and 
dividing  the  sum  by  the  new  ship  weight.  The  equation  for  KG  is: 

KGnew  =  (KGold  x  Wold  +  KG  add  x  Wadd)  ■*■  Wnew 
The  Simulink  representation  is  shown  in  Figure  5. 

F.  VERTICAL  CENTER  OF  BUOYANCY  (KB) 

The  ship's  new  vertical  center  of  buoyancy  (KB)  is  determined  by  multiplying  the 
original  KB  by  the  original  displacement,  adding  the  product  of  the  added  weight  and  the 
original  draft  plus  one  half  of  the  parallel  sinkage,  then  dividing  the  product  by  the  new 
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Figure  4:  SIMULINK  iMoment  Terminology 
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Figure  5:  SIMULINK  KG  Terminology 
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weight.  In  equation  form: 

KB  NEW  =  (KBORIC  X  WORIG  +  KB  ADD  X  WaDd)  +  WnEW 

This  is  based  on  the  assumption  that  the  hull  is  effectively  wall-sided,  and  that  the  vertical 
center  of  buoyancy  of  the  newly  submerged  portion  is  located  at  one  half  the  added  draft. 
Both  of  these  assumptions  are  consistent  with  the  small  angle  assumption  that  was 
mentioned  in  the  previous  chapter.  If  the  small  angle  assumption  is  no  longer  true,  then 
the  true  position  of  the  center  of  buoyancy  will  have  to  be  calculated.  This  can  be  done  by 
employing  numerical  integrations  and  utilizing  the  provided  table  of  offsets,  or  by 
interfacing  the  program  with  a  standard  external  ship  hydrostatics  calculations  program, 
such  as  SHCP  [Ref.  2:  pp.  24-25].  Such  an  interface  is  relatively  straight  forward  since 
the  computational  environment  of  Matlab/Simulink  already  allows  incorporation  of 
external  ASCII  files.  Figure  6  shows  the  Simulink  representation  of  the  KB  equation. 

G.    HEIGHT  OF  THE  METACENTER  ABOVE  THE  CENTER  OF  BUOYANCY 
(BM) 

The  height  of  the  metacenter  above  the  center  of  buoyancy  (BM)  is  determined  by 

the  second  moment  of  the  waterplane  (about  the  waterplane  centerline  or  longitudinal 

center  of  floatation)  divided  by  the  hull  displacement.  For  the  added  weight  method,  a 

free  surface  correction  is  made  by  subtracting  the  second  moment  of  the  flooded 

compartment  area  from  the  waterplane  moment.  Typically,  free  surface  corrections  are 

done  for  liquids  when  calculating  the  final  metacentric  height.  In  such  a  case,  the  free 

surface  correction  appears  as  a  reduction  in  the  metacentric  height  or,  equivalently,  an 

increase  in  the  virtual  center  of  gravity.  The  same  effect  can  be  observed  by  modifying 
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the  second  moment  of  the  waterplane  area  in  the  calculation  of  the  metacentric  radius, 
BM.  This  is  because  the  added  weight  in  this  case  is  seawater.  With  this  modification  in 
BM,  the  added  weight  method  produces  identical  results  with  the  lost  buoyancy 
calculations  (reference  a  naval  architecture  book).  The  BM  Simulink  terminology  is 
shown  in  Figure  7. 

H.  RIGHTING  ARM  (GM) 

The  righting  arm  (GM)  is  determined  by  GM  =  KB  +  BM  -  KG.  This  equation  is 
graphically  represented  by  the  triple  summing  point  shown  in  Figure  8. 

I.      HEEL  AND  TRIM  ANGLE 

The  angle  of  heel  is  determined  by 

sin  0  =  Volume  x  Ztrans  -s-  (Displacement  x  GM ) 
as  shown  in  Figure  9.  The  total  trim  is  determined  by  dividing  the  total  trimming  moment 
by  the  moment  to  change  trim  (MCT).  MCT  is  equal  to  the  total  weight  times  the 
longitudinal  GM,  divided  by  the  ship  length  at  the  waterline.  The  tangent  of  the  trim 
angle  is  equal  to  the  total  trim  divided  by  the  ship's  length  at  the  waterline.  If  the  large 
angle  assumption  needs  to  be  utilized,  the  moment  to  change  trim  will  have  to  be 
imported  as  an  external  array.  This  is  not  expected  to  cause  significant  difficulties  since 
MCT  is  readily  available  from  the  hydrostatic  diagrams  of  a  given  ship.  The  Simulink 
trim  angle  terminology  is  shown  in  Figure  10. 
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J.  ADDITIONAL  SUBMERGENCE 

The  additional  submergence  of  the  hole  due  to  list  is  determined  by  multiplying  the 
sine  of  the  heel  angle  by  the  transverse  distance  from  the  transverse  center  of  floatation  to 
the  hole.  Similarly,  the  increased  hole  depth  due  to  trim  is  obtained  by  multiplying  the 
sine  of  the  trim  angle  by  the  longitudinal  distance  to  the  center  of  floatation.  The  new 
hole  depth  is  the  sum  of  the  original  hole  depth  and  the  depth  increases  due  to  heel  and 
trim.  In  algebraic  form,  these  relations  are  shown  by 

hlNEw  =  h2oLD  +  (Ztrans  . hok  x  sin  <p )  +  (Zlong  , hole  x  sin(  trim  )) 

and  in  graphical  form  they  are  represented  by  the  diagram  of  Figure  1 1 . 


K.  PROGRAM  USAGE 

To  use  Pflood.m,  the  user  needs  only  select  "simulation"  from  the  SIMULINK 
menu,  then  select  "start".  Modifications  to  the  simulation  parameters  or  the  program  are 
easily  accomplished  by  "double-clicking"  on  the  symbol  of  interest.  Additionally, 
SIMULINK  allows  the  user  to  easily  add  output  displays  in  order  to  observe  the  behavior 
of  any  desired  variable  during  the  simulation. 
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IV.  RESULTS 

Results  for  a  test  case  consisting  of  a  compartment  with  a  4  square  foot  hull 
penetration  (compartment  1)  and  a  4  square  inch  watertight  bulkhead  penetration 
(compartment  2).  Pumps  were  activated  at  10  minutes  (600  seconds)  which  completely 
dewatered  compartment  1   in  order  to  test  the  ability  of  the  program  to  account  for 
backflow  through  holes.  As  shown  in  figure  12,  the  four  square  foot  hull  penetration 
completely  floods  compartment  1  essentially  instantly.    Activating  the  pumps 
dewaters  compartment  1 ,  instantly  removing  the  water  entering  through  the  hull 
penetration.  Similar  results  could  have  been  achieved  by  changing  the  characteristics  of 
the  hull  penetration  to  simulate  hull  repairs. 

Figure  1 3  illustrates  how  the  program  accurately  predicts  backflow  into 
compartment  1  from  compartment  2,  which  had  no  pumps  assigned.  A  saturation  function 
ensures  that  backflow  from  compartment  2  will  stop  when  compartment  2  water  height 
drops  to  the  level  of  the  bulkhead  penetration. 

The  predicted  trim  response  of  the  ship  is  shown  in  Figure  14.  The  ship  initially 
trims  as  a  result  of  the  rapid  flooding  of  compartment  1 ,  then  gradually  increases  its  trim 
angle  as  compartment  2  floods.  The  instant  dewatering  of  compartment  1  then  causes  the 
trim  angle  to  decrease,  followed  by  further  trim  angle  reduction  as  compartment  2  is 
dewatered.  Without  altering  the  pump  or  hole  configuration,  the  ship  will  assume  a  steady 
trim  angle  corresponding  to  compartment  2  flooded  up  to  the  height  of  the  bulkhead 
penetration. 
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Figure  12:  Height  of  Water  in  Compartment  1 
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Figure  13:  Height  of  Water  in  Compartment  2 
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Figure  14:  Ship  Trim  Response 
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V.  DISCUSSION,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A.  DISCUSSION 

This  thesis  has  successfully  developed  the  foundation  for  a  modular  and 
expandable  progressive  flooding  computer  model.  The  computer  model  accurately 
predicts  the  quasi-static  ship  response  to  progressive  flooding  and  is  easily  modified  by 
the  user  to  explore  ship  response  to  any  desired  damage  situations. 

B.  CONCLUSIONS 

Based  on  the  above  results,  the  following  conclusions  are  drawn: 

1 .  Pflood.m  is  an  effective  design  tool  for  modeling  and  studying  the  effects  of 
progressive  flooding. 

2.  Pflood.m  is  easily  modified  to  model  virtually  any  damage  scenario. 

3.  Pflood.m  can  be  expanded  to  include  different  ship  progressive  flooding  and 
equilibrium  models. 

C.  RECOMMENDATIONS 

The  following  items  are  recommended  to  further  expand  Pflood.m  and  enhance 
its  utility  as  a  design  tool: 

1.  Incorporate  a  hull  characteristics  look-up  table  into  the  Pflood.m  program 
to  expand  the  ship  response  modeling  ability  beyond  "small  angles"  of 
heel  and  trim  and  non-rectangular  compartments. 

2.  Explore  the  effect  of  pump  configuration  on  dewatering  effectiveness  and 
ship  response. 
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22.922 
23.074 
23 .212 
.587  . 
1.965 
3.209 


20.066 

24.0675 

28.069 

29.591 

31.113 

031 

.2815 

.532 


4.329  .782 
5.37  1.032 
7.259  1.5325 
8.943  2.033 
11.835  3.034 

079 

833 
.197 
.271 

105 

258 

0  12 


4.035 

5.038 

6.041 

7.045 

8.049 

10.055 

.061 


10. 

10. 

10 

10. 

10. 

10 

10 

10 

10, 

10 

10 

10 


14 

15 

17 

IS 

19 

20 

21 

21.301  13.238 

21.566  14.415 

22.0513  17.2415 

22.442  20.068 

22.827  24.069 

23.177  28.07 

23.253  29.3665 

23.318  30.663 

.587  .031 

1.9484  .2815 

3.189  .532 

4.329  .782 

5.39  1.032 

7.322  1.5325 

9.062  2.033 


10.0 

10.0 

10 

10 

10 

10 

10 

10 

10 


12.085 

14.498 

16.407 

17.896 

19.051 

19.944 

21.167 

21.915 

22'.  197 

22.422 

22.761 

23 

23 


034 
035 
0375 
04 
044 
8.048 
10.0545 
12.061 
13.239 
14.417 
17.2435 
0  20.07 
232  24.07 


23.392  28.072 
23.4195  29.0725 
23.442  30.073 
784  .587  .031 
784  1.8596  .2815 
784  3.034  .532 
10.784  4.13  .782 
10.784  5.15  1.032 
10.784  7.0292  1.5325 
10.784  8.733  2.033 
10.784  11.725  3.034 
10.784  14.164  4.035 
10.784  16.126  5.037 
10.784  17.704  6.039 
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10.784  13.969  7  043 
10.784  19.953  8  047 
10.784  21.267  10  0535 
10.784  22.075  12  06 
10.784  22.381  13.2385 
10.784  22.612  14  &1 
10.784  22.943  17  2*435 
10.784  23.142  20  07 
10.784  23.362  24  071 
10.784  23.482  28.072 
10.784  23.494  28  9075 
10.784  23.502  29  743 
11.961  .588  .031 
H-961  1.674  .2815 
11.961  2.69  .532 
11.961  3.647  .7825 
11-961  4.548  1.033 
11.961  6.222  1.533 
11.961  7.764  2.033 
11.961  10.515  3.034 
11.961  12.917  4.035 
11-961  15.02  5.0365 
11.961  16.768  6.038 
11-961  18.183  7.041 
11-961  19.286  8.044 
11-961  20.846  10  05 
11.961  21.839  12  05 
11.961  22.219  13  237 
11-961  22.49  14.416 
11-961  22.899  17  243 
H-961  23.122  20.07 
'  H-961  23.372  24.07 
H-961  23.502  28.07 
11.961  23.509  28.7325 
11.961  23.512  29.393 
12.941  .589  .031 
12.941  1.435  .282 
12.941  2.258  .533 
12.941  3.043  .783 
12.941  3.8  1.033 
12.941  5.229  1.533 
12.941  6.575  2.034 
12.941  9.042  3.034 
12.941  11.294  4  035 
12-941  13.369  5.036 
12.941  15.211  6.037 
12-941  16.806  7.0395 
12.941  18.118  8.042 
12.941  20.04  10.0485 
12.941  21.258  12  0^5 
12.941  21.732  13.2345 
12.941  22.081  14  414 
12.941  22.574  17  24 
12.941  22.882  20.069 
12.941  23.222  24.07 
12.941  23.462  28.071 
12.941  23.488  28.657 
12.941  23.512  29  243 
14.02  .59  1.001 
14.02  2.1554  1.517 
14.02  3.66  2.033 
14.02  6.326  3.034 
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14. 

02  8 

.563  4.035 

14. 

02  10.76  5.0355 

14. 

02  1 

.2.739  6.035 

14. 

02  1 

.4.602  7.0375 

14. 

02  16.244  8.039 

14. 

02  18.732  10.045 

14. 

02  20.266  12.05 

14. 

02  20.87  13.23 

14. 

02  21.316  14.41 

14. 

02  21.981  17.239 

14. 

02  22.432  20.067 

14. 

02  22.962  24.068 

14. 

02  23.382  28.069 

14. 

02  23.437  28.641 

14. 

02  23.492  29.213 

15. 

098 

.59  2.941 

15. 

098 

2.318  3.487 

15. 

098 

4.091  4.033 

15. 

098 

6.998  5.0335 

15. 

098 

9.546  6.034 

15. 

098 

11.783  7.035 

15. 

098 

13.83  8.036 

15. 

098 

17.1296  10.042 

15. 

098 

18.954  12.048 

15. 

098 

19.658  13.228 

15. 

098 

20.189  14.409 

15. 

098 

21.081  17.2365 

15. 

098 

21.731  20.064 

15. 

,098 

22.551  24.065 

15. 

,098 

23.261  28.066 

15. 

,098 

23.393  28.6595 

15. 

.098 

23.493  29.25 

16, 

.078 

.588  5.067 

16. 

.078 

4.447  6.05 

16. 

.078 

7.762  7.033 

16. 

.078 

12.144  8.535 

16 

.078 

15.153  10.037 

16 

.078 

17.441  12.045 

16 

.078 

18.677  14.053 

16 

.078 

19.138  15.056 

16 

.078 

19.524  16.059 

16 

.078 

20.769  20.06 

16 

.078 

22.004  24.06 

16 

.078 

22.7355  26.722 

16 

.078 

23.433  29.383 

16 

.863 

.587  6.801 

16 

.863 

3.559  7.417 

16 

.863 

6.27  8.032 

16 

.863 

12.98  10.036 

16 

.863 

15.967  12.04 

16 

.863 

16.83  13.04 

16 

.863 

17.42  14.05 

16 

.863 

17.955  15.05 

16 

.863 

18.419  16.056 

16 

.863 

19.969  28.058 

16 

.863 

21.34  24.06 

16 

.863 

22.236  26.821 

16 

.863 

23.074  29.583 

17 

.647 

.587  8.541 

17 

.647 

10.03  10.23 

17 

.647 

14.02  11.92 

17 

.647 

15.292  12.98 
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17 

.647 

16.08  14.046 

17 

.647 

16.74  15.05 

17 

.647 

17.298  16.054 

17 

.647 

18.958  20.0565 

17 

.647 

20.39  24.059 

17 

.647 

21.375  26.951 

17 

.647 

22.309  29.843 

18 

.039 

.586  9.406 

18 

.039 

7.97  10.467 

18 

.039 

12.014  11.528 

18 

.039 

13.203  12.11 

18 

.039 

14.027  12.708 

18 

.039 

14.931  13.553 

18 

.039 

15.649  14.399 

18 

.039 

17.331  17.229 

18 

.039 

18.399  20.058 

18 

.039 

19.805  24.059 

18 

.039 

21.15  28.06 

18 

.039 

21.464  29.027 

18 

.039 

21.774   29.993 

18 

.431 

.585  10.276 

18 

.431 

4.648  10.704 

18 

.431 

8.007  11.132 

18 

.431 

10.22  11.583 

18 

.431 

11.64  12.034 

18 

.431 

13.443  13.039 

18 

.431 

14.63  14.044 

18 

.431 

15.474  15.049 

18 

.431 

16.087  16.053 

18, 

.431 

17.798  20.057 

18. 

.431 

19.165  24.06 

18, 

.431 

20.171  27.1065 

18. 

,431 

21.148  30.153 

18, 

,824 

-585  11.146 

18. 

,824 

5.835  11.589 

18. 

824 

9.241  12.032 

18. 

824 

11.024  12.453 

18. 

824 

12.018  12.875 

18. 

824 

13.29  13.635 

18. 

824 

14.191  14.395 

18. 

824 

15.54  16.2255 

18. 

824 

16.414  18.056 

18. 

824 

17.498  21.058 

18. 

824 

18.47  24.06 

18. 

824 

19.45  27.194 

18. 

824 

20.419  30.328 

19. 

118 

-584  11.797 

19. 

118 

5.067  12.08 

19. 

118 

8.005  12.36 

19. 

118 

10.783  12.863 

19. 

118 

12.018  13.365 

19. 

118 

12.914  13.879 

19. 

118 

13.565  14.392 

19. 

118 

14.338  15.2215 

19. 

118 

14.906  16.051 

19. 

118 

15.903  18.0555 

19. 

118 

16.64  20.06 

19. 

118 

17.91  24.0605 

19. 

118 

19.12  28.061 

19. 

118 

19.48  29.282 

19. 

118 

19.834  30.503 

19. 

608 

.584  12.881 
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19. 

508 

5.626    13.076 

19. 

608 

8.004    13.271 

19. 

608 

10.317    13.5625 

19. 

608 

11.33    13.854 

19. 

608 

11.944    14.122 

19. 

608 

12.402    14.39 

19. 

608 

13.355    15.269 

19. 

608 

14.037    16.147 

19. 

608 

15.747    20.104 

19. 

608 

16.936    24.061 

19. 

608 

17.853    27.39 

19. 

608 

18.713    30.719 

19. 

804 

.584    13.311 

19. 

804 

4.002    13.406 

19. 

804 

6.002    13.501 

19. 

804 

8.3687    13.671 

19. 

804 

10.006    13.841 

19. 

804 

11.4077    14.149 

19. 

804 

12.024    14.457 

19. 

804 

13.25    15.63 

19. 

804 

14.038    16.804 

19. 

804 

15.479    20.433 

19. 

804 

16.507    24.062 

19. 

804 

17.412    27.443 

19. 

804 

18.213    30.823 

20. 

0     .584    14.38 

20. 

0    4. 

4    14.561 

20. 

0    6. 

005    14.741 

20. 

.0    8. 

273    15.0815 

20. 

.0    1C 

1.01    15.422 

20. 

.  0    1] 

..75    15.82 

20. 

.o  i: 

S.339    16.217 

20. 

.  0    12 

1.73    16.916 

20 

.  0    14 

1.038    17.615 

20 

.0    14.534    18.8375 

20 

.0    14.95    20.06 

20 

.0    16.1    24.062 

20 

.0    17.08    28.064 

20 

.0    17.4    29.4955 

20 

.0    17.703    30.927 

20 

.049 

.584    15.542 

20 

.049 

2.815    15.618 

20 

.049 

4.003    15.694 

20 

.049 

5.975    15.896 

20 

.049 

7.402    16.097 

20 

.049 

9.225    16.421 

20 

.049 

10.764    16.746 

20 

.049 

12.442    17.158 

20 

.049 

13.821    17.57 

20 

.049 

14.58    19.278 

20 

.049 

15.126    20.987 

20 

.049 

15.619    22.631 

20 

.049 

16.049    24.276 

20 

.049 

16.465    25.939 

20 

.049 

16.852    27.602 

20 

.049 

17.215    29.276 

20 

.049 

17.563    30.95 

20 

.093 

.584    16.586 

20 

.093 

3.194    16.677 

20 

.093 

4.375    16.768 

20 

.093 

6.534    17.014 

20 

.093 

8.138    17.26 
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20. 

093 

10.138  17.654 

20. 

093 

11.85  18.049 

20. 

093 

13.102  18.407 

20. 

093 

14.29  18.7658 

20. 

093 

14.824  20.251 

20. 

093 

15.268  21.736 

20. 

093 

15.772  23.571 

20. 

093 

16.237  25.405 

20. 

093 

16.667  27.258 

20. 

093 

17.062  29.11 

20. 

093 

17.252  30.04 

20. 

093 

17.437  30.971 

20. 

098 

.584  16.7 

20. 

098 

4.488  16.881 

20. 

098 

6.005  17.061 

20. 

098 

8.304  17.402 

20. 

098 

10.009  17.742 

20. 

098 

11.04  17.972 

20. 

098 

12.01  18.202 

20. 

098 

13.264  18.523 

20. 

098 

14.322  18.849 

20. 

098 

14.543  19.455 

20. 

098 

14.75  20.06 

20. 

098 

15.35  22.062 

20. 

098 

15.89  24.063 

20. 

098 

16.598  27.017 

20. 

098 

17.23  29.97 

20 

103 

.584  16.819 

20 

103 

2.903  16.896 

20 

103 

4.055  16.97 

20 

103 

6.03  17.18 

20 

103 

7.504  17.388 

20 

103 

9.364  17.724 

20 

103 

10.912  18.06 

20 

103 

12.757  18.507 

20 

103 

14.267  18.954 

20 

103 

14.86  20.584 

20 

.103 

15.37  22.214 

20 

.103 

15.836  23.89 

20 

.103 

16.26  25.565 

20 

.103 

16.655  27.262 

20 

.103 

17.008  28.959 

20 

.196 

.584  19.032 

20 

.196 

'2.629  19.1 

20 

.196 

3.779  19.168 

20 

.196 

5.634  19.343 

20 

.196 

6.956  19.519 

20 

.196 

8.664  19.806 

20 

.196 

10.1  20.09 

20 

.196 

11.796  20.478 

20 

.196 

13.204  20.862 

20 

.196 

13.976  21.074 

20 

.196 

14.746  21.285 

20 

.196 

15.05  22.114 

20 

.196 

14.34  22.943 

20 

.196 

15.78  24.397 

20 

.196 

16.17  25.851 

20 

.196 

16.569  27.4 

20 

.196 

16.947  28.951 

20 

.294 

.584  21.361 

20 

.294 

4.413  21.536 

20 

.294 

6.005  21.711 
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20.294 

20.294 

20.294 

20.294 

20.294 

20.294 

20.392 

20.392 

20.392 

20.392 

20.392 

20.392 

20.392 

20.392 

20.392 

20.392 

20.392 

20.49 

20.49 

20.49 

20.49 

20.49 

20.49 

20.49 

20.49 

20.49 

20.49 

20.49 


8.225  22.049 
10.009  22.386 


13.064 
15.414 
16.149 
16.884 
.588  23 


23 
23 
26 
28 
.6' 


1275 
869 
406 
943 


4.155  23.851 
5.975  24.031 
8.307  24.3715 
10.009  24.712 


12.188 
14.012 
14.953 
15.762 
16.284 
16.805 

.534  26. 

3.117  26 


25.2275 
25.743 
26.0365 
26.33 
27.637 
28.943 
Oil 
086 


4.004  26.161 
7.627  26.599 
10.009  27.037 


12.125 
13  .929 
15.086 
16.131 
16.506 
16.88 


27 
28 
28 

28 


535 
033 
4075 
782 


28.8625 
28.943 
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APPENDIX  B.  ARVOL.M 


arvol .m 

by  Russell  wright 

calculates  station  areas  and  center  o '  bouvancy 

for  T  =16.9  ft 


load  hulll .dat 
a=zeros (41,3) ; 
zz=16.9 
for  i=l:41 

beg=23* (i-l)+l; 

endd=23*i; 

a(i, l)=hulll(beg)  ; 

z=hulll (beg: endd,  3  )  ; 

if  z (1)  >  zz 

a(i,3)=0,-a(i,2)=l; 
else 
%   y=hulll(beg:endd,2) 

zr)=hulll  (beg:  endd,  3  ) 
z=zp(l) : (16.9-zp(l) ) /22:16.9, 
yy= hulll (beg: endd,  2) 
y=SDline (zD,yy, z)  ; 
h=z"(2)-z(l)  ; 
hhh(i)=h; 

four  =  0.0;fr=0; 
two=0 .0; tO=0; 

for  j=  1:11 

four=four+y (2* j ) ,- 

fr=fr+z(2*j)*y(2*j) ; 
end 

for  k=l:10 

two-  two+y ( 2*k+l) ; 

to=to+y(2*k+l) *z(2*k+l) ; 
end 

%  calculate  station  area  using  Simpson's  rule: 

a(i,2)=2* (y(l)+4*four+2*two+y(23) ) *h/3; 

%  calculate  station  centroid  (z  bar)  using  Simpson's  rule: 

a(i,3)=2* (y(l)*z(l)+4*fr  +  2*to+y(23) 'z (23 ) ) *h/ ( 3*a ( i , 2 )  )  ; 

end 

end 

aa=zeros  (101 ,  3  )  ,- 

hh=( (a(37,l)-a(3,l)  )  /100)  ; 

aa  ( 

aa( 

aa  ( 


, l)  =  (a(3,l) :hh:a(37,l) )  ' 
,2)=spline(a(3:37,l)  ,a(3 
,3)=spline(a(3:37,  1)  ,a(3 


37,2) ,aa( : ,1) ) 
37,3) ,aa( : ,1) ) 


f ore=0 ; too=0 ; f re=0 ; tw=0 ; f o  =  0 ; tu=0  , 
for  i=l:50 

f ore=fore+aa(2*i, 2) ; 

fre=f re+aa(2*i, 2) *aa (2*i,  3  )  ; 

fo  =  fo+aa(2*i,2) *aa(2*i,  1)  ; 
end 

for  i=l:49 
too=too+aa(2*i+l , 2) ; 
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tw=tw+aa(2*i+l,2)*aa(2*-!+]    3) 
tu=tu+aa(2*i+l,2) *aa(2*i+~'i) 
end 


vol=(hh/3)*(aa{1.2)+4*fore+2*too+aa(101   2))- 
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APPENDIX  C.  DRAFTS.M 


%% 

%  draf ts2 .m 

%  by  russell  wright 

%  calculates  wp  points  for  a  given  draft 

% 

load  hulll .dat 
ind=0; 

yincl=zeros (41,3) ; 
zz=16.9; 

zincl=zz 
zinc2=zz 
for  i=l:41 
i 

beg=23* (i-l)+l; 

endd=23*i; 

wpla(i, 1) = hulll (beg, 1) ; 

z=hulll (beg: endd, 3 ) ; 

y= hulll (beg: endd, 2 )  ; 

wpla ( i , 2 ) =interpl (z ,y, zz, ' spline' ) ; 
if  z ( 1)  >  zz 

wpla(i,2)=0; 

yincl (i,2)=0; 

yincl (i, 3) =0; 

wpla(  i,  1)  =0  ,• 
end 

yincl (i, 1) =hulll (beg,  1)  ; 

yincl ( i , 2 ) =interpl (z,y, zincl, ' spline' ) ; 

if  yincl (i,2)  <0 .0 

yincl (i, 2)  =0.0; 
end 
yincl ( i , 3 ) =zz; 

if  wpla(i, 2) ==0 
yincl (i, 2) =0; 
yincl (i, 3 ) =0; 
end 

end 

y  inc  1  =y  inc  1(3:37,:),- 

xx=yincl ( : , 1 ) ; 

x^yincl (1,1)  :  (yincl (35, 1) -yincl (1,  1)  ) / 34 : yincl ( 3  5 , 1)  ; 

yy=yincl  (  :  ,  2 )  ,- 

y=spline (xx,yy , x) ; 

h=X ( 2 ) -x ( 1 ) ; 
hhh(i)=h; 

four  =  0.0;fr=0; 
two=0  .  0  ;  to=0  ,- 

for  j=  1:17 

f our=four+y (2' j  )  ; 
.  fr=fr+x(2*j)*y(2*j) ; 
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end 

for  k=l:16 

two=  two+y (2*k+l) ; 

to=to+y (2*k+l) *x(2*k+l)  ; 
end 

%  calculate  wp  1/2  area  using  Simpson's  rule: 
a2=2*(y(l) +4*four+2*two+y (35) ) *h/3; 

%  calculate  wp  centroid  (x  bar)  using  Simpson's  rule: 
xbar=2*  (y(l)  *x(  1)  +4*fr+2-to+y (35)  *x(35)  )  *h/  (3*a2) 

%  calculate  long,  moment  o'  inertia  about  xbar 


f  r=0  ; 
to=0; 

for  j=  1:17 

fr=fr+(x(2*j)~2)*y(2*j)  ; 

end 

for  k=l:16 

to=to+y(2*k+l) * (x(2*k+l) "2) ,- 


io=2*(y(l) *x(l) *x(l)+4*fr+2*to+y(35)*(x(35) ~2) )  *h/3; 
isubl=io-a2*xbar~2 
%  calculate  trans  moment  o'  inertia 


fr=0; 
to=0; 

for  j=  1:17 

fr=fr+(y(2*j)*3)  ; 

end 

for  k=l:16 

to=to+(y(2*k+l)~3) ; 
end 

isubt=2* ( (y(l)~3)+4*fr+2*to+(y(35)~3) ) *h/9; 
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